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Abstract 

This paper describes state-of-the-art digital filtering techniques that are part of GEANA, an automatic data analysis 
software used for the GERDA experiment. The discussed filters include a novel, nonlinear correction method for ballistic 
deficits, which is combined with one of three shaping filters: a pseudo-Gaussian, a modified trapezoidal, or a modified cusp 
hlter. The performance of the filters is demonstrated with a 762 g Broad Energy Germanium (BEGe) detector, produced 
by Canberra, that measures 7 -ray lines from radioactive sources in an energy range between 59.5 and 2614.5 keV. At 
1332.5 keV, together with the ballistic deficit correction method, all filters produce a comparable energy resolution of 
~ 1.61 keV FWHM. This value is superior to those measured by the manufacturer and those found in publications 
with detectors of a similar design and mass. At 59.5 keV, the modified cusp filter without a ballistic deficit correction 
produced the best result, with an energy resolution of 0.46 keV. It is observed that the loss in resolution by using a 
constant shaping time over the entire energy range is small when using the ballistic deficit correction method. 

Keywords: Digital Hltering techniques, ballistic deficit correction, high purity germanium detectors (HPGe), radiation 
detection, 7 -ray spectroscopy, GERDA 


1. Introduction 

The use of High Purity Germanium (HPGe) detec¬ 
tors is reliant on the application of effective signal filter¬ 
ing methods to properly estimate the amount of energy 
deposited by radiation. Notably, the evaluation of the 
data quality for experiments with multiple detectors re¬ 
quires fast and accurate digital signal processing routines. 
GEANA (’’GErmanium ANAlysis”) is a novel digital sig¬ 
nal processing software created to address this need in the 
scope of the GERmanium Detector Array (GERDA) ex¬ 
periment [T]. It is capable of monitoring the steady flow 
of data collected by the experiment over several years. 

The energy deposited by radiation in the germanium 
detectors is reconstructed by processing the digitized sig¬ 
nals with digital filters. Digital filters compensate for 
effects, such as the decay of the signal induced by the 
RG-feedback of the charge sensitive preamplifier, the fi¬ 
nite charge collection time of the signal and frequency 
components dominated by noise. This publication aims 
to describe the filters used to counter these effects within 
GEANA. 

Shaping filters are used to improve the signal-to-noise 
ratio. The pseudo-Gaussian filter, the most common digi¬ 
tal shaping filter, brings the input pulse close to a Gaussian 
distribution. This shaping filter is fast in execution, simple 
to implement and its performance is sufficient for many ap¬ 
plications. However, for particular situations other filters 
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perform more optimally. The current signal of a typical 
HPGe detector, unfortunately, does not have the shape of 
a delta function. Nevertheless, if it is assumed that a delta 
function is an input signal to a circuit that only adds cur¬ 
rent and voltage noise, the infinite cusp filter PJ p. 245], 
has the best performance. However, the infinite nature 
of the cusp filter conflicts with the finite length of digital 
signals. 

The time profile of the current signal (pulse shape), 
that is induced into the read out electronics varies with 
the number and location of energy deposition sites in a 
HPGe detector. Neither the cusp nor the pseudo-Gaussian 
filter are well suited to deal with variations in the pulse 
shape. These hlters underestimate the amount of energy 
deposited in the detector for pulses with long rise times 
relative to pulses with a short rise time. This effect is 
called ballistic deficit and can be partially compensated 
for by choosing a sufficiently large shaping time or the 
trapezoidal filter [3]. 

In the past, two additional approaches were proposed 
to correct ballistic dehcits with analog modules. The first 
approach calculated a semi-empirical compensation which 
depended on the input signal rise time and then was added 
to the final output signal |3]. The second approach used 
a gated integrator after prefiltering the signal to integrate 
over the charge collection time [S]. 

2. Methods 

The different steps of the pulse height (energy) recon¬ 
struction in GEANA are illustrated in Fig. [T] and discussed 
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Figure 1: The different steps of the pulse height reconstruction 

applied to a signal with a pileup. The three initial steps illustrate 
the pole-zero cancellation. It is followed by a differentiation and the 
shaping filter (pseudo-Gaussian with a shaping time of 6 ps). During 
the differentiation the analysis range is reduced to a window around 
the pulse that is wide enough to not considerably affect the shaping 
filter’s performance. 



Figure 2: The effect of the multi-site event cancellation (MSEC) 
method on three different current pulses (normalized). All the con¬ 
tributions around the pulse are added up to a single sample (the size 
of that samples exceeds the y-range up to roughly 1.) 


with d = exp(A/r), where A is the sampling width and r 
the decay time of the preamplifier. This inverse high-pass 
filter can be applied to the signal to compensate the decay 
introduced by the RC-feedback circuit. However, due to 
an unknown offset of the baseline, the output signal will 
have a nonzero slope (see second trace of Fig. This 
slope can be removed by a simple linear interpolation. The 
interpolation is performed by dividing the baseline (the 
recorded samples in the range before the charge pulse) into 
two segments and calculating their averages. From the two 
averages the offset and slope of the baseline is extracted 
and subtracted from the entire pulse. As shown in Fig.[2 
this procedure corrects also distortions from the decaying 
tail of proceeding events (pileups). 


in the following sections. First, the decay introduced by 
the RC-feedback circuit of the preamplifier in the recorded 
signal (charge pulse) is compensated with a pole-zero can¬ 
cellation. This signal is differentiated (current pulse) and 
can be corrected for ballistic deficits with a nonlinear, ir¬ 
reversible filter, the multi-site event cancellation (MSEC). 
Finally the shaping filter is applied. The shaping filters are 
based on a combination of three linear 0 p. 87] filters: the 
high-pass, low-pass, and moving average filters; their re¬ 
cursive implementations are discussed on p. 277 and p. 322 
of Ref. [5] . 

2.1. Pole-zero cancellation 

A high-pass filter has the same impact on a step-like 
pulse as the RC-feedback circuit does. The recursive form 
of the high-pass filter can be inverted (to calculate the n-th 
output sample from the input samples Xi): 

, . d r. 1 ~ ^ \ 

yn = (g-Xn+>x^)-m, g=- --, m = 2——, (I) 

^ I - a 1 -b d 


2.2. Ballistic deficit correction 

The multi-site event cancellation (MSEC) method has 
been developed to correct for variability in the pulse shape. 
It is particularly tailored to correct for events with energy 
depositions in multiple locations (multi-site event, MSEs) 
that present particularly strong variations in the charge 
collection time. It is similar to the method that uses a 
gated integrator, however, it does not prefilter the signal. 
Instead, after pole-zero correction and differentiation, the 
signal is set to zero in a window that stretches over the 
duration of the charge collection process. The subtracted 
contributions are then added to the position, where the 
pulse previously reached its maximum. The width of the 
window to which the filter is applied is found indepen¬ 
dently for each individual signal. This is done by search¬ 
ing for the first sample to the left and to the right of the 
maximal value of the current signal that falls below a cer¬ 
tain threshold. The threshold level is extracted from the 
baseline of the signal. The altered signal is then fed into 
the shaping filter. The respective procedure is applied to a 
variety of pulses in Fig. to show the effect of the MSEC 
method. 
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2.3. Shaping Filters 

The pseudo-Gaussian filter is composed of multiple con¬ 
secutive moving average filters of the same width [ 6 l p. 281]. 
The impulse response of the first moving average filter 
takes the shape of a rectangular function, the second mov¬ 
ing average filter changes this into a triangular function. 
Each further moving average filter transforms the impulse 
response closer to a Gaussian shape. In this work the 
pseudo-Gaussian filter is composed of four consecutive mov¬ 
ing average filters. 

The standard trapezoidal filter is constructed with two 
moving average filters of different width. The difference in 
width of the two moving average filters defines the length 
of the flat-top of the filters output signal. The pulse height 
ideally is reconstructed at a constant position towards the 
end of the flat-top [7]. However, such a location is not 
easily found for a pulse with a variable time profile. This 
ambiguity is bypassed by applying two additional moving 
average filters (of half the size of the flat-top) to change 
the flat-top to a rounded-top and form a clear, distin¬ 
guished maximum. The width of the rounded-top is de¬ 
fined through the width of the initial flat-top portion (the 
difference in the width of the two initial moving average 
filters). 

The standard cusp filter can be constructed through 
a bidirectional [SI p. 330] low-pass filter. In this publica¬ 
tion an additional moving average filter flattens the peak 
of the traditional impulse response. This creates a round- 
top shape, similar to the one described in Ref. | 8 ]. This is 
better suited to the finite charge collection time of HPGe 
signals. The width of the rounded-top is defined by the 
width of the additional moving average filter. The cusp 
filter optimally extends to infinity, however, in a discrete 
implementation its width must be limited. Thus, the im¬ 
plemented version is applied in the range of ±5 times the 
shaping time around the charge pulse. 

The altered versions of the two filters will be referred 
to as rounded-top trapezoidal and rounded-top cusp filter. 



Figure 3: The 1332.5 keV 7 -ray line of ®®Co (pseudo-Gaussian filter, 
shaping time of 6 ps, ballistic deficits corrected) and the different 
components of the peak shape model. The uncertainty of the FWHM 
considers only the uncertainty of the fitting procedure. 


duced by incomplete charge collection and ballistic deficits: 
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The variables cr, Xg, (3, the number of events in the peak 
Sg and the tail St, and the two amplitudes A, B are free 
parameters of the model. After completing the fitting pro¬ 
cedure, the energy resolution is calculated from the model 
as the full-width at half-maximum (FWHM) value of the 
the Gaussian peak and the decaying tail, i.e. gix) tix). 
Fig. I^shows the ®°Co line at 1332.5 keV together with the 
model described in this section. 


2.Jt. Evaluation of the energy resolution 

Statistical fluctuations in the number of produced charge 
carriers and electronic noise broaden a 7 -ray line’s peak 
into a Gaussian shape. The width of the broadening is 
measured through the energy resolution, which is a funda¬ 
mental characteristic of a radiation detector. 

There are different ways to define the energy resolution 
and it is crucial to clearly define the respective procedure. 
Within this publication the energy resolution is evaluated 
by fitting a function (model) to the shape of the peak in the 
pulse height spectrum by using the Levenberg-Marquardt 
technique [5]. The model m consists of a Gaussian peak 
gix), a step-like shelf s(a:), and a Hypermet function tix). 
The latter is added to account for a possible tail, which 
decays exponentially below the peak’s centroid and is pro¬ 


3. Setup 

To examine the functionality and performance of the 
methods used in GEANA, data were collected with a stan¬ 
dard Canberra Broad Energy Germanium (BEGe) detec¬ 
tor m- The detector is made from p-type germanium, 
has a mass of 762 g and a cylindrical form with a radius of 
3.6 cm and a height of 3.5 cm. A small read-out electrode 
is embedded in one base of the detector and most of the 
remaining surface is covered by a high voltage electrode. 
The BEGe detector is mounted in a 7500SL vertical dip¬ 
stick vacuum cryostat, which contains a 2002CSL charge 
sensitive preamplifier (~ 50 ps decay time). The output 
of the preamplifier is amplified to fit the dynamic range 
of the analog to digital converter. The acquisition system 
consists of a Struck SIS3301 VME flash analog-to-digital 
converter (FADG) [H]. The FADG module accommodates 
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Figure 4: A scatter plot of the reconstructed energy and the rise 
time distribution around the 1332.5 keV 7 -ray line of ®®Co with the 
multi-site event cancellation (MSEC) enabled and disabled (pseudo- 
Gaussian filter, shaping time of 6 |xs). 


up to eight input channels and is equipped with an anti¬ 
aliasing bandwidth filter. The signals are sampled in 40 ns 
steps and quantized with 16 bit resolution. The length of 
the signal is defined to be 164 ps and the rising edge of the 
pulse is in the center of the sampled time interval. 

Two measurements were taken to test the filters, one 
with a dynamic energy range up to 1750 keV using a va¬ 
riety of uncollimated radioactive sources (^"^^Am, 
i34Cs, i37Cs, 6 OC 0 , 228te) and one with an energy 

range up to 3000 keV with a uncollimated ^^®Th source. 


4. Results 

The ®°Co 7 -ray line at 1332.5 keV is a standard 7 -ray 
line in 7 -ray spectroscopy and is used for most of the dis¬ 
cussion. Where useful, lines at low energies (59.5 keV of 
^^^Am, 121.8 keV of ^^^Eu) and at high energies (2614.5 keV 
of ^^®Th) are also studied to show differences in a broad en¬ 
ergy range. Unless otherwise stated, the standard configu¬ 
ration is the pseudo-Gaussian filter with a shaping time of 
6 ps and MSEC enabled. The following discussion is valid 
for BEGe-like detectors. Nevertheless, the MSEC method 
can also be used to correct for variations in the pulse shape 
found in other detector types. 

4.I. Performance of the MSEC method 

Fig.g] shows a scatter plot of the energy reconstructed 
with a pseudo-Gaussian filter of 6 ps shaping time and the 
rise time for the region around the 1332.5 keV 7 -ray line of 
®°Co. In the left panel ballistic deficits are not corrected. 
The shaping filter therefore reconstructs a lower energy for 
events with a long rise time. In the right panel ballistic 
deficits are corrected with the MSEC method, resulting 
in the energy peak being reconstructed for any rise time 
value in a much narrower energy interval. 



Figure 5: The energy resolution (FWHM) of the 1332.5 keV 7 -ray 
line of ®''Co obtained with the pseudo-Gaussian filter for different 
shaping times. The uncertainties are not drawn; their values are less 
than 0.5% and thus not visible. 
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Figure 6 : The best energy resolutions (FWHM) found individually 
for each 7 -ray line with the pseudo-Gaussian filter (top) and the cor¬ 
responding optimum shaping time (middle). The difference in energy 
resolution of an analysis with a constant shaping time optimized for 
the 1332.5 keV 7 -ray line and the best energy resolution found for an 
individual 7 -ray line is represented in the bottom. 
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In addition, the shaping time of the pseudo-Gaussian 
hlter can be varied. Fig. [^illustrates, how the energy res¬ 
olution depends on both, the shaping time and the MSEC 
method at 1332.5 keV. The resolution is dominated by 
voltage noise at short shaping times and by current noise 
at long shaping times. The optimum shaping time is found 
at the position where the energy resolution is minimal and 
the two contributions are equal (the noise corner). A simi¬ 
lar behavior can be observed for any 7 -ray line in the pulse 
height spectrum. 

Effects from ballistic deficits further broaden the en¬ 
ergy resolution and shift the optimum shaping time and 
the minimal reachable energy resolution to higher values. 
The MSEC method corrects this effect and improves the 
optimum shaping time from 12 ps to 6 ps at 1332.5 keV 
(see Fig. [^. The energy resolution at that energy addi¬ 
tionally improves by roughly 4% from (1.669 ± 0.009) keV 
to (1.612 ± 0.006) keV. A short shaping time allows for a 
measurement at higher count rates. Next to improving the 
energy resolution, the MSEC correction also reduces the 
optimum shaping time and increases the maximum event 
rate that can be processed. 

4-S. Energy dependence of the MSEC method 

The top panel of Fig. [^ presents the curve that re¬ 
lates the energy resolution of the most prominent peaks 
to the energy. At low energies the resolution is primar¬ 
ily dominated by electronic noise. With increased energy 
the contribution from statistical fluctuations in the num¬ 
ber of charge carriers increases and dominates the energy 
resolution. MSEs from Compton scattering are the most 
important reason for variations in the pulse shape in BEGe 
detectors and thus the most important cause of ballistic 
deficits. In germanium, Compton scattering is surpassed 
by photoelectric absorption as the dominant interaction 
channel below 100 — 200 keV. With the disappearance of 
Compton scattering at low energies the effects caused by 
ballistic deficits also disappear. 

Above ^ 400 keV the MSEC method improves the en¬ 
ergy resolution on average by roughly 50 eV. The improve¬ 
ment is not energy dependent. Below ~ 400 keV the im¬ 
provement quickly decreases and below 200 keV the MSEC 
method actually deteriorates the energy resolution. The 
reason for this is that the summing operation of the MSEC 
method allows more noise to pass through the shaping fil¬ 
ter in respect to the situation where the MSEC method is 
not applied. 

If the ballistic deficits are not corrected then the opti¬ 
mum shaping time is energy dependent. Ballistic deficits 
grow with increasing energy which in turns increases the 
optimum shaping time to compensate the ballistic deficits 
(see Fig. [^ at these energies. This effect is illustrated in 
the central panel of Fig. [^ The effect is reduced consid¬ 
erably if the MSEC is enabled, but not entirely removed. 
This suggest that the MSEC method is capable of remov¬ 
ing an important fraction of the effects caused by ballistic 
deficits but not all of them. 


In a typical analysis the shaping time optimization is 
not performed for each peak individually, instead, a single 
shaping time is chosen for analyzing the entire energy spec¬ 
trum. Most of the peaks therefore are not reconstructed 
with the optimum shaping time and the best energy res¬ 
olution is not always obtained. An analysis with a con¬ 
stant shaping time has been performed for both situations, 
MSEC disabled and enabled. The shaping time that was 
used to reconstruct the spectrum was selected by finding 
the optimum value at the 1332.5 keV 7 -ray line of ®°Co 
(12 ps for the MSEC disable and 6 ps for the MSEC en¬ 
abled). The energy resolutions from a spectrum obtained 
with only a single shaping time are then compared to the 
energy resolutions that were found when the shaping time 
was optimized (by finding the noise corner) at each peak 
individually. The bottom panel of Fig. [^ presents the 
difference of the energy resolution obtained with a sin¬ 
gle shaping time and the value found at the noise cor¬ 
ner of each line. The maximal deviation when the MSEC 
method is disabled is roughly 150 eV. When the MSEC 
method is enabled the loss in energy resolution is below 
40 eV. This is also valid if the peaks from ^^®Th, which 
are above 1500 keV, are included. Accordingly, with the 
pseudo-Gaussian filter and the MSEC method enabled it is 
possible to analyze the entire energy range from 59.5 keV 
up to 2614.5 keV with a single shaping time without sig¬ 
nificant losses in the energy resolution. Therefore, the 
MSEC method is particularly important for an analysis 
that only uses a single shaping time over a large range of 
energies. This is the case for most applications in 7 -ray 
spectroscopy. 


4-3. Trapezoidal and cusp filter 

The MSEC method can also be combined with other 
filters such as the rounded-top cusp and the rounded-top 


trapezoidal filter introduced in Sec. 2.3 


At 1332.5 keV, with MSEC enabled the two filters have 
a similar performance to the pseudo-Gaussian filter for all 
rounded-top widths (see the upper plot of Fig. MSEC 
on). On average, the tested filters found a resolution of 
1.611 keV with a standard deviation of 0.004 keV. Thus, 
the differences between the filters are within the uncertain¬ 
ties of the respective fits. Also the same optimum shaping 
time is found within a narrow range of values. Where bal¬ 
listic dehcits dominate, the choice of a specihc filter does 
not matter much as long as MSEC is enabled. 

The picture is different if ballistic deficits are not cor¬ 
rected and the real sensitivity of the filters to these ef¬ 
fects appears (see the upper plot of Fig. [7[ MSEC off). 
At 1332.5 keV with MSEC disabled filters with a wide 
rounded-top perform better as more of the ballistic deficits 
are removed. Moreover, all filters have a tendency towards 
higher shaping times as this also helps to remove ballistic 
deficits. The rounded-top cusp hlter is least suited to be 
used in the presence of strong ballistic dehcit effects and 
even with a large rounded-top width it is not capable of 
matching the performance of the other two hlters. The 
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Figure 7: The best energy resolution and the optimum shaping 
time of the 1332.5 keV (top) and the 59.5 keV (bottom) 7 -ray lines 
for the rounded-top cusp and rounded-top trapezoidal filter with var¬ 
ious rounded-top widths (see Sec. |2.3[ |. For comparison, the pseudo- 
Gaussian filter is shown. 

pseudo-Gaussian filter performs slightly better, however, 
only the rounded-top trapezoidal filter reaches a perfor¬ 
mance comparable to the situation with MSEC enabled. 
This requires a large rounded-top width but allows for a 
much lower shaping time than observed with the other hl- 
ters. 

At 59.5 keV the situation is quite different (see bottom 
plot of Fig. as ballistic deficits are not present. Similar 
to the pseudo-Gaussian hlter, the energy resolution found 
with the two other filters is worse if the MSEC method 
is enabled. When the MSEC method is disabled, at these 
energies the rounded-top cusp hlter reaches its full noise 
reduction potential and performs better than the other 
two hlters. Both Hlters with an adjustable rounded-top 
perform best when the rounded-top width is small. 

The hlter’s behavior is tightly linked with both, the 
presence of ballistic dehcits and their performance in re¬ 
ducing noise. When ballistic dehcits are present it is vital 


to correct for these effects. Best suited to do this is the 
MSEC method, but a rounded-top trapezoidal hlter can 
also be used. In the absence of ballistic dehcits (i.e. at low 
energies) a cusp hlter (with a small rounded-top width) 
reduces the noise more efficiently than the other two hl¬ 
ters and thus performs better. However, it is important to 
note, that because of its slow falling form, the cusp hlter 
requires a much larger signal range and is not suitable to 
high count rate measurements. 

4.4- Comparison with other measurements 

Some of the evaluated energy resolutions can be com¬ 
pared to the values provided by the manufacturer of the 
BEGe detector. The manufacturer quoted an energy reso¬ 
lution of 1.77 keV at 1332.5 keV, which is higher than the 
averaged value of (1.611 ± 0.004) keV found in the previ¬ 
ous section. The value measured here is also better than 
the resolution of ~ 1.63 keV found in Ref. [13] and the 
best value of (1.65 ± 0.01) keV found in Ref. [14]. Both 
publications used a very similar detector for their stud¬ 
ies. The lowest FWHM value measured in this study was 
(0.582 ± 0.001) keV for the ^®^Eu y-ray line at 121.8 keV. 
This value can be compared to the energy resolution of 
0.66 keV measured by the manufacturer at an energy of 
122.1 keV with a ^"^00 source. These comparisons show 
that the implemented filter methods perform well. 

The best energy resolution obtained for the 7 - 

ray line at 59.5 keV is (0.463 ± 0.001) keV and the 
7 -ray line at 2614.5 keV is (2.315 ± 0.008) keV. The indi¬ 
cated uncertainties are the standard deviation of the fit 
parameters. 

5. Conclusion 

Important aspects of digital signal processing techniques 
used in radiation detectors, such as different shaping fil¬ 
ters and a novel method to correct ballistic deficits have 
been examined for a Canberra Broad Energy Germanium 
detector (BEGe). With the multi-site event cancellation 
(MSEC) method enabled all filters lead to comparable re¬ 
sults at 1332.5 keV, the differences are within the uncer¬ 
tainties. The trapezoidal filter with a large round-top has 
the best performance at 1332.5 keV if ballistic deficits are 
not corrected. At energies below ^ 200 keV a ballistic 
deficit correction is not required, instead, it is more bene¬ 
ficial to use a cusp like filter at these energies to find the 
best signal-to-noise ratio. 

The optimum shaping time for a specific 7 -ray line has 
been found to be energy dependent. But even an analysis 
with a single shaping time over a large energy range did 
deviate by less than 40 eV (pseudo-Gaussian filter) from 
the optimal value, provided the MSEC is applied. Energy 
peaks above a few hundred keV are dominated by Comp¬ 
ton scattered events and a correction for ballistic deficits 
is necessary to avoid an even larger deterioration of energy 
resolution. 
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These filters are implemented into GEANA. The soft¬ 
ware can filter up to SO’OOO events per second (on a single 
Xeon E5-2660 2.2GHz GPU), if the data can be supplied 
with an I/O of over 300MB/s. Furthermore, the software 
also performs a fully automated energy calibration and fil¬ 
ter optimization. GEANA has already successfully been 
used as an independent cross check in the Phase I data 
release m of GERDA. No discrepancies have been found 
in the measured energy range, which reaches up to 7MeV. 
Furthermore, it was extensively used for the pulse shape 
analysis of semi-coaxial detectors in that data release [16] . 
In the second data collection phase of the experiment, 
GEANA is being adopted to analyze data gathered from 
photomultiplier tubes and silicon photomultipliers that are 
part of the new liquid argon instrumentation m- 
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